Majorana liquids: the complete fractionalization of the electron 
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We describe ground states of correlated electron systems in which the electron fractionalizes 
into separate quasiparticles which carry its spin and its charge, and into real Majorana fermions 
which carry its Fermi statistics. Such parent states provide a unified theory of previously studied 
fractionalized states: their descendants include insulating and conducting states with neutral spin 
5* = 1/2 fermionic spinous, and states with spinless fermionic charge carriers. We illustrate these 
ideas on the honeycomb lattice, with field theories of such states and their phase transitions. 



o 

> 

O 

O 
o 



I 

o 
o 



> 

in 

^' 

o 
o 



The study of two-dimensional quantum antiferromag- 
nets has proved to be fertile ground for finding many- 
electron states whose excitations do not carry all the 
quantum numbers of the electron This phenomenon 
is often referred to as "spin-charge" separation. It leads 
to some of the most non-trivial examples of quantum 
entanglement at long scales, and is crucial for the under- 
standing of a variety of correlated electron materials, and 
for designing topological quantum computers. 

Upon fractionalizing the electron into its spin and 
charge, one is faced with the decision of locating its Fermi 
statistics. In early theories of gapped spin liquid states 
of insulating, frustrated, antiferromagnets, distinct phys- 
ical motivations led to two main pictures: 
[i) A picture of projecting out doubly-occupied and va- 
cant sites from a free electron Slater determinant pro- 
duced an attachment of Fermi statistics to spin, leading 
to spin liquids with neutral, spin S = 1/2 excitations 
('spinous') which are fermions. 

{ii) A picture of quantum-'disordering' magnetically or- 
dered states produced an attachment of Fermi statistics 
to charge, leading to ,3] spin liquids with bosonic spinous. 
This dichotomy persists when we consider insulating spin 
liquids with gapless excitations, and also to metallic or 
superconducting fractionalized states. Thus in the first 
picture we have 'algebraic spin liquids' (ASL) with Dirac 
oints or Fermi surfaces of neutral fermionic spinous 0- 
while in the second picture we have gapless bosonic 
spinous and/or gapless, spinless fermionic charge carriers 
[9l4ll ■ including 'algebraic charge liquids' (ACL). 

This paper will provide a unification of these seem- 
ingly divergent pictures. We will argue that the states 
above descend from a common parent 'Majorana liquid', 
in which there is complete fractionalization of the elec- 
tron into components carrying its spin, charge, and Fermi 
statistics. No choice is made in the parent liquid regard- 
ing attachment of Fermi statistics, which is instead car- 
ried by real Majorana fermions which carry neither spin 
nor charge. While we will introduce our approach gener- 
ally, a specific realization of the theory will be provided 
on the honeycomb lattice. We will also relate our results 
to a recent numerical study of the Hubbard model on the 
honeycomb lattice [l2| . 

We motivate our approach by starting from the sec- 



ond picture Here, we imagine there is some pre- 

ferred local magnetic order, and transform the electron 
(ccr, <7 =tii) into a rotating reference frame determined 
by the orientation of the local magnetic order [l3| : 
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Here the spinful bosons Zo- define a SU(2) rotation ma- 
trix TZz determined by the magnetic order, and /i, /2 are 
spinless fermions which carry the charge of the electron. 
It is also useful to note the various global and gauge sym- 
metries associated with Eq. ([T|) . The global spin rotation 
SU(2)spin acts as a left multiplication on TZz, while the 
global charge U(l)chargc is carried by /i_2- In addition 
there is a local SU(2)s^g gauge invariance: this acts as a 
right multiplication of TZ^ and the doublet (/i, /2) trans- 
forms as its fundamental. 

Now let us turn to the first picture [3, 01. We can 
view this as a transformation into a rotating reference 
frame in the pseudo-spin space of particle-hole transfor- 
mations of the electron 14]. For the Hubbard model on 
bi-partite lattices at half-filling, there is in fact a global 
SU(2)pseudospin which contains U(l)chargo as a subgroup; 
however our approach only assumes U(l)charse,S]/_mmetry 
in general. Paralleling Eq. ([T]) we now have [; 
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Independent of the existence of a global SU(2)psoudospin 
symmetry, this parameterization introduces a local 
SU(2)c.g gauge invariance: this acts as a right multi- 
plication of T^b and the doublet {fi,fl) transforms as 
its fundamental. This SU(2)c,g gauge invariance played 
a crucial role in the classification of various spin liquid 
states with fermionic spinous (l5. 16|. 

There is a natural and simple unification of Eqs. ([1]) 
and ([2]). We write the complex fermions in terms of 2 
sets of real fermions Ca and Xa, with a = 1 . . . 4, by c-|- = 
Ci + *C2 and cj, = ^3 + id and similarly between /i_2 and 
Xa- Then we have 
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where 7?. is a real S0(4) matrix. This shows that a com- 
bination of Eqs. ([1]) and ([2]) enjoys a S0(4)g gauge in- 
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variance; it will become evident below that S0(4)g ~ 
SU(2)3,g (g) SU(2)c,g. As before, the S0(4)g gauge in- 
variance acts as a right multiplication of 7?., x is a fun- 
damental of S0(4)g, and global symmetries act as left 
multiplication of TZ. For an explicit form of the global 
symmetries, let us represent the 4x4 real matrices as 
tensor products of 2 sets of 2 x 2 Pauli matrices: s" act- 
ing on the t, i space (a — x, y, z), and acting on the 
Re[c], Im[c] space. Then the global SU(2)spin symmetry 
is generated by 

s"" = s^py , sy = sy , s*^ = (4) 

The global pseudo-spin transformations are generated by 

= syp'- , Ty = s^p^ , = py- (5) 

here generates U(l)chargc, while T^^y generate the re- 
maining pseudo-spin transformations which need not be 
a symmetry of the Hamiltonian. The matrices in Eqs. ^ 
and ([5]) are the 6 generators of SO (4); note also that the 
3°^ commute with the , realizing the factorization into 
SU(2) ® SU(2). To complete our formulation, we can ex- 
press TZ in terms of the complex bosons in TZz and TZb- 
Writing = (l)Q — i(f)^, = —(j)2~i4>\, and bi — (ftQ + icf)'^, 
62 = ^*02 + '/'i where (/)^ and 0^ are real scalars, we have 

z, ^ (\,i + i^\T^ + i(\qry + i^%T\ (6) 

Eqs. dni) and ^ contain our general statement of elec- 
tron fractionalization: the electron ^ decomposes into the 
bosonic fields Z^ and Z^ which carry its global SU(2)spin 
and U(l)charge quantum numbers respectively, and into 
the Majorana fermion x carrying the Fermi statistics. 
The resulting theory has a S0(4)g = SU(2)s,g ® SU(2)c,g 
gauge invariance: Z^ and x carry SU(2)s_g charges, and 
Zc and X carry SU(2)c,g charges. 

Different patterns of breaking the S0(4)g gauge invari- 
ance and global symmetries lead to a plethora of possible 
phases, and we present a broad classification: 

(I) Phases with conventional excitations are obtained 
when both and Zs are condensed. In such phases, we 
can always choose a S0(4)g gauge in which Zc = Za = 1, 
and then it becomes clear from Eq. ([3]) that the fermion 
X has just the same quantum numbers as the electron 
C. By condensing various fermion bilinears (just as 
in conventional Hartree-Fock/BCS theory), we can ob- 
tain Fermi liquids, semi-metals, antiferromagnets, va- 
lence bond solids (VBS), superconductors, or quantum 
spin Hall states [rzl-list. 

(II) When Zg is condensed, we can use the SU(2)s,g gauge 
invariance to set Zs = \. Then Eq. ((3|) reduces to Eq. dZl), 
and we therefore reproduce the phases of Refs. on 
the square lattice, and those of Ref. @ on the honeycomb 
lattice, including the ASLs. In these phases, the fixing of 



the SU(2)3^g gauge transfers the global SU(2)spin quan- 
tum numbers from Z^ to Xj while the U(l)chargc quantum 
numbers remain on Z^.. Thus these phases have neutral 
fcrmionic spinons, and bosonic charge carriers. 

(III) A complementary situation is realized when Z^ is 
condensed. Now we can set Z^ = \, Eq. ^ reduces to 
Eq. ([1]), global U(l)chargo quantum numbers are trans- 
ferred from Zc to Xj and so we obtain neutral bosonic 
spinons and fermionic charge carriers. We reproduce 
phases of Refs. [13, Ei on the square lattice, and Ref. [13] 
on the honeycomb lattice, including the ACLs. 

(IV) When neither Zf. and Zg are condensed, we can ob- 
tain phases in which the Z^, Z^^ and x are separate ele- 
mentary excitations, carrying the charge, spin, and Fermi 
statistics of the electron respectively. These are the Ma- 
jorana liquids of this paper. These elementary excitations 
all carry gauge quantum numbers, and so the stability 
of such phases requires that gauge forces not be confin- 
ing: we will describe specific examples of deconfinement 
mechanisms below. 

The remainder of the paper applies our general the- 
ory to the half-filled, extended Hubbard model on the 
honeycomb lattice. For weak interactions, we have a 
conventional semi-metal with electronic excitations at 2 
Dirac points in the Brillouin zone: this is as realized in 
graphene. For str ong interactions, there is convincing nu- 
merical evidence [l2| for an insulator with coUinear, two- 
sublattice antiferromagnetism (Neel order) . The simplest 
possibility is that there is a direct transition between the 
these two category I phases [20| . However, recent numer- 
ical studies [12] indicate there may be an intermediate 
phase. 

We begin our analysis by describing the free electron 
spectrum in the semi-metal phase, associated with the 
Hamiltonian 

i/o = -<5]cU<, (7) 

where i are sites on the honeycomb lattice, and (ij) refers 
to nearest neighbors. We take the low energy limit of 
in the standard manner, obtaining two valleys of two- 
component Dirac fermions. Explicitly, we expand the 
electron at two Dirac valleys by di,2 = e*'^^'^'''c (where 
(5i,2 = i(x'^) wavevectors of the valleys), and 

introduce Pauli matrices r" and pF' which act on the 
sublattice and valley spaces respectively. Then, after 
introducing real Majorana fermions as the real and 
imaginary parts of e^T"^ e'^'^ (di, 1x^(^2)*, we obtain the 
continuum Lagrangian 

8 

^O^^CaTM^MCa. (8) 
0=1 

Here /i is a 2-1-1 dimensional spacetime index, and the 
Dirac 7 matrices are (70,71772) = T^,T^,r^. In contrast 
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to the single-site Majorana fermion in Eq. ([3]), here the 
fermion field has additional components associated with 
the sublattice and valley spaces. The sublattice index is 
equivalent to the spacetime Dirac index and has been left 
implicit, and from now on the 'flavor' index a — 1 ... 8 
accounts for the spin, pseudospin, and valley indices; thus 
for each a, (^a is now a 2-component Majorana spinor. We 
can now decompose ^ as in Eq. 1^, with Xa having the 
same spinor structure as Ca, while TZ is as in Eq. ([6]). 

We obtain our parent algebraic Majorana liquid (AML) 
when interactions beyond those in Hq leave both Zs and 
Zc un-condensed and realized as gapped quanta which 
carry spin and charge respectively. The AML is in 
our category IV above, and at energies below the spin 
and charge gaps, it has gapless, relativistic, Majorana 
fermions Ca coupled to emergent SU(2) gauge fields 
and A"^ associated with the SU(2)s^g SU(2)c,g gauge 
invar iance: 
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X7m (^m 



lAl^T^) X (9) 



The stability of the AML requires that the gapless x 
fermions suppress monopoles, and so that the SU(2)s.g ® 
SU(2)c,g gauge forces are not confining 0, [2l|. Such 
monopole suppression happens for a sufficiently large 
number of gapless fermion flavors, and it is not known 
if the 8 real fermion flavors here are sufficient. Assum- 
ing deconfinement, the AML has a gap to all spin and 
charge excitations, and has gapless Majorana fermions 
which carry only energy. 

Whether or not the AML is stable, we can use it to 
describe a very large number of descendant phases. The 
rest of the paper will note some interesting or physically 
relevant examples. 

First, we can expect that the SU(2)s,g®SU(2)c,g gauge 
forces lead to the analog of chiral symmetry breaking, 
and the simplest possibility is the appearance of a 0(8) 
invariant XX condensate. Such a condensate leads to a 
fermion mass gap, and breaks time-reversal symmetry, 
leading to a chiral Majorana liquid , also in category IV. 
The fermions generate Chern-Simons terms for the gauge 
fields, and this leads to deconfinement 22| for the gapped 
Zc, Zs, and x excitations. There is non-zero spin chiral- 

ity S\ ■ {S2 X S3) -I , and also nonzero electrical currents 

on the lattice. However, the physical current is sup- 
pressed from XX- 



(CO - {xx){zhzlz,.jZc,,)«: 



(10) 



i and j are two next nearest neighbor sites, therefore 
the expectation value {Zl ^Z\ ^Zs.jZc.j) is expected to be 
small when Zs and Zc are both gapped. 

There are also a large number of possible Higgs phases. 
One interesting example is the Higgs condensate of the 
vector if": 



If H" had involved bilinears of the original electron 
C, its condensate would break spin rotation invariance 
and lead to a quantum spin Hall phase In the 

present situation, the H" in Eq. (fTT|) does not carry 
any global quantum numbers, and its Higgs condensate 
does not break any global symmetries. The resulting 
phase is in fact a Z2 Majorana liquid, and is in category 
IV, as we now argue. The if" condensate breaks the 
SU(2)s_gC3>SU(2)c,g gauge invariance to U(l)s_g(g)U(l)c,g; 
if we choose H" oc (0, 0, 1), then the U(l)'s are generated 
by and T^. Thus the low energy theory of this phase 



X7m i^tJ. 



mxp^S^Xi 

(12) 

where the fermion mass term is induced by the H°' con- 
densate. We can further integrate out the massive x 
fermions, and then using the analog of the arguments 
in Ref. 18|, we find that the physics is controlled by a 



mutual Chern-Simons term 
2i 



r — 



2^ tf^P c^^i^v^s.p- 



(13) 



As discussed elsewhere [2J], with such a term, the gauge 
forces are quenched, and the matter fields carry only Z2 
gauge charges. The Z2 gauge field endows mutual statis- 
tics between excitations with and charges, for in- 
stance between the spin- and charge-carrying bosons, Zg 



and Zc, as in Ref. |23l |. 



The last two category IV phases above, the chiral and 
Z2 Majorana liquids, are attractive candidates for the in- 
termediate state in Ref. [l^ between the semi-metal and 
the Neel insulator. They have gapped Zc, Zg, and x exci- 
tations, and we have demonstrated that the charge, spin, 
and Fermi statistics they carry remain deconfined. These 
phases treat the charge and spin excitations at an equal 
footing (unlike the proposal in Ref. [3!), and so they are 
appropriate for the vicinity of the metal-insulator transi- 
tion. The chiral Majorana liquid has weak spontaneous 
spin chirality and electrical currents, which have not been 
detected so far. The Z2 Majorana liquid has no broken 
symmetry, and so remains compatible with all existing 
computations. 

Let us now turn to category II. As noted earlier, such 
phases have (Zg) ^ 0, and on the honeycomb lattice 
our theory reduces to that of Hermele [6|]. He found 
an SU(2)c,g ASL insulator: at low energies, this is de- 
scribed by the theory of the AML in Eq. but with 
the SU(2)s,g gauge fields A"^ = because of the Zg 
condensate. As with the AML, this is stable only for suf- 
ficiently large number of fermion flavors. The x fermions 
now carry spin (explained in (II)), and so the spin exci- 



H''=x P^'S'^ x- 



(11) 



tations are gapless. The intermediate phase of Ref. [12 1 
has a spin gap in addition to the charge gap, and this is 
not compatible with an ASL. 

Finally, we turn to category HI, where we have (Zc) ^ 
0. These were discussed in Ref. ll| for the square lattice, 



4 



(A) 

{Zs) + 
(7V)^0 
Neel 

(B) 

{Z,) = 

VBS 



(C) 

(^.) ^ 
(7V) = 
semi-metal 



(D) 

[Z,) = 



(7V) = 
SU(2),,, 



: 

ACL 



FIG. 1: Some phases of the half-filled honeycomb lattice. All 
phases above have the condensate (Zc) 7^ 0, and so can be de- 
scribed by Eq. Phase D is category III and is described 
by Eq. @ but with = 0: it has a spin gap, and gap- 

less, spinless critical fermions at the Dirac points carrying 
the electronic charge. Phases A and C are conventional and 
so in category I: C has gapless electron-like excitations at the 
2 Dirac points, which are gapped in A. Phase B is described 
by Eq. H15|) . which undergoes monopole- induced confinement 
to VBS order. 



As after Eq. p2|) . we can integrate out the massive 
fermions, but now find only a Maxwell term for the A], ^ 
gauge field. Thus the low energy theory of phase B has 
a gapless, relativistic U(l)s^g photon. In the absence of 
gapless matter, it is known that the monopoles in such 
a gauge field condense, and lead to long-range order de- 
termined by the quantum numbers of the monopole [25| . 
We will describe the computation of monopole quantum 
numbers elsewhere, showing that a kekule type valence 
bond solid (VBS) order develops. The same conclusion 
is reached by approaching phase B from phase A [25| . 
Thus, while phase B started out in category III, it ulti- 
mately becomes category I VBS state. 

This paper has unified two previously divergent ap- 
proaches to the study of the electron fractionalization: 
those with fermionic [1, El, [iBl versus bosonic (ol-flU 
spinous. We have applied the theory to the honeycomb 
lattice and predicted new phases of possible relevance to 
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and we generalize the discussion here to the honeycomb 
lattice. As in category II above, we begin with the AML 
in Eq. ©, and now derive a SU(2)s_g ACL obtained by 
setting the SU(2)c,g gauge fields A"^ — 0. This phase 
has a spin gap, but the fermions x '^O'^ form gapless 
excitations which carry charge (as explained in (III)), 
again incompatible with the phase of Ref. [l2|- As for 
the AML, the ACL can lead to a 'chiral charge liquid' 
with a XX condensate. 

Other phases which descend from the SU(2)s^g ACL 
are shown in Fig. [T] The simplest of these is the semi- 
metal phase C: this is the category I phase with both 
and Zs condensates, as noted earlier. For the remaining 
phases in Fig. [T] we have to consider a Higgs condensate 
of the field 

TV" = X m''^" X- (14) 

This choice is motivated by the fact that if we replace x 
by the electron ^, then Eq. (|14p is the conventional Neel 
order; this replacement is permitted when both Z^ and 
Zs are condensed. This is the case in phase A, which is 
then a category I phase with Neel order. The transition 
between phases C and A involves the order parameter 
and the gapless Dirac electrons, but no gau ge f ields; it is 
in the class of field theories studied in Ref. [20|. 

Now we consider the remaining phase B in Fig.[T] By a 
gauge transformation, we orient the iV" condensate along 
(0,0,1); such a condensate breaks the SU(2)s_g gauge 
invariance down to U(l)s,c,, leaving only the A\^^ gauge 
field active. The neutral Zs spinons are gapped, and so 
the low energy theory of phase B (following Eq. ([12])) is 

-Cb = X7m (^m - ^K,^S') X + m XM^^'X- (15) 
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